Stationary flows in quantum dissipative closed circuits as a 
challenge to thermodynamics 



V. Capek and J. Bok 
Institute of Physics of Charles University, 
Faculty of Mathematics and Physics, 
Ke Karlovu 5, 121 16 Prague 2, Czech Republic 
(Tel. (**-420-2)2191-1330, Fax (**-420-2)2492-2797, E-mail capek@karlov.mff.cuni.cz) 

(Dated: April 24, 2002) 

Abstract 

An experimentally inspired model is constructed and rigorously solved from the Hamiltonian 
level where a dc circular spontaneous flow exists in absence of a magnetic field, irrespective of 
presence of dissipation causing otherwise proper relaxation. The flow causes a spontaneous uni- 
directional transfer of heat from one bath to another one, even against temperature step. This 
is what is explicitly forbidden by the Clausius form of the Second law of thermodynamics. The 
unidirectionality of the flow is caused by that of spontaneous processes known to bear this property 
since their introduction by Einstein. The model slightly improves the previous one (Capek &: Shee- 
han 2002), describes a realistic plasma system for which experimental results violating the second 
law were announced, and the result obtained fully supports the experimental conclusions (Sheehan 
1995). Analytical proof of the violation is supported by numerical results. All mathematical details 
are exposed, two fully independent types of mathematical arguments behind starting equations are 
invoked, and no approximations that could be made responsible for the striking conclusions are 
used. It shows how the physics beyond the Second law is still little understood. 

PACS numbers: 05.30.-d, 05.70.-a, 44.90. +c 
Keywords: 2"'' Law Challenge 
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I. INTRODUCTION 



Challenges to the second law of thermodynamics |T], are almost as old as the law itself, 
dating back at least to the 1870's with Maxwell's celebrated demon 0, ^, ^ 0. Most of 
them have been resolved under close scrutiny ^, ^ but some persisted. Anyway, strong 
belief in old authorities and natural human tendency to organize things and facts into 
closed logical units and complexes (scientific disciplines etc.) caused that almost nobody 
doubts about validity of standard thermodynamics, in particular the Second law, in at least 
the macroworld |TU|. Absolutistic statements like '...No exception to the second law of 
thermodynamics has ever been found - not even a tiny one...' [TU| often appear and the 



second law is almost universally believed to be unquestionable. In our opinion, the situation 
is (in view of the fast developing situation perhaps at least still) not so clear. The reason 
is that experiments questioning the second law have been reported since 1995, have since 
been subject to a public discussion but remain so far unquestioned Jll], |12[- Some of these 
experiments have since been even reproduced (compare |0 with [Q). 



In theory, arguments independent of these experiments appeared since 1997 |]T5[ saying 
that in quantum systems with strong or at least intermediate coupling to its surroundings 
(identical or connected with usual thermodynamic baths) with mutual strong correlations 
(entanglement), the standard statistical thermodynamics could be violated. This is in par- 
ticular, but not only, the case of the second law of thermodynamics. One must keep in mind 
that from first principles (microscopic Hamiltonian dynamics), derivations of the second law 
are declared to exist just in classical (in the sense of non-quantum) physics []T^ or, as in 
standard textbooks, involving assumption of a weak system-bath coupling only. ^ It should 
be stressed already here that the classical physics is, according to the Bohr correspondence 
principle, an infinite temperature limit of the (more general) quantum physics. This, inter 
alia, implies that its application to finite temperatures as in standard thermodynamics is at 
least open to discussion. The above lack of general derivation of the second law beyond, in 
particular, the classical regime could also correspond to the fact that so far reported and 
seemingly classical paradoxes connected with the second law (see, e.g., [|I^) usually involve 



^ Gibbs canonical form of the system density matrix compatible with the standard thermodynamics is 
correct to just the zeroth order in the system-bath coupling. This is an important fact to be realized 
already here as the effect reported below is of higher order in this coupling. 
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sufficiently intense processes that are inherently of the quantum character. 

In 1999, it was realized that a long lasting call in chemistry for inclusion of self- 
organizational tendencies into theory of particle-transfer chemical reactions is, from the 
microscopic point of view, nothing but a call for inclusion of such mechanisms that can turn 
any (from the thermodynamical point of view) passive bath into an active one, opening thus 



door to violations of standard thermodynamic principles |T3. Recent review of theoretical 
models and state-of-art in theory could be found in, e.g., |]19|, |20|. So far, two main groups of 
theoretical models of purely quantum open systems violating the second law existed: Those 
with quantum reaction channels opening or closing in accordance with the instantaneous 
state of the system (reminding of the Maxwell demon closing and opening gate in a 
wall separating two compartments with a classical gas), and those where a specific type 



of interference of different quantum reaction channels exists [19|. In particular this type 
of systems is relevant as the contradiction with the second law treated in [|19| is not only 
mathematically well justified but can be given even a very simple physical interpretation 
based on otherwise experimentally well established facts: Exciton diffusion bearing energy 
and going (as always diffusion does) in the direction of decreasing exciton concentration. 



Recently, as a theoretical response to another positively tested experimental system 11 



a next model of still another type has been suggested ^T|] that also allows rigorous solution 



(exact within a scaling theory) fully conffiming experimental doubts about universal validity 
of the second law. Detailed discussion and solution of a modified (and nearer to reality) 
version of the model is the subject of the present paper. Here, however, we are already able 
to support the conclusions by both scaling and non-scaling arguments. What is perhaps 
universal for all the models challenging thermodynamics is that the system in question must 
be, during its activity, outside the canonical state p2|^ . Mechanisms how to achieve and 



maintain that might be, of course, different. One should understand, however, that this 
condition is perhaps necessary but by far not sufficient. For other paradoxical systems that 
could be also classified as above, showing how physics beyond the second law is still little 



^ A comment might be relevant in this connection that outside equihbrium, there is no reason why the 
thermodynamic entropy (characterizing individual macroscopic body) and the Shannon - von Neumann 
entropy (characterizing ensemble, i.e. bearing another physical information) should be identical. For 
validity of the second law in Nature, the former entropy is relevant while only the latter is usually 
involved in analytical 'proofs' of the second law starting from first principles. 
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understood see, e.g., O, E3 . 



II. THEORETICAL MODEL 

Physicists are more conservative than mathematicians. For the latter, one unquestioned 
proof is standardly enough to accept a new statement, theorem etc. Physicist on the other 
hand require usually more arguments, especially when old truths or dogmas are challenged. 
That is why we should like to stress already here that: 

• The model treated here is by no means unique. It is just a next one in a long series 
of models behaving, according to a strict microscopic theory, in a way incompatible 
with standard thermodynamics. Among those, it is exceptional just in the sense of 
modelling an experimental system tested positively, from the point of view of violability 
of the second law, already a long time ago. 

• Except for model assumptions, we shall use no approximate steps in our reasoning in 
the main text here. This is to be stressed in order to avoid misunderstandings. 



Mathematics we are going to use in the main text is that one by Davies ^ which 
forms, for one specific choice of the scaling parameter, mathematical basis of the weak- 
coupling kinetic relaxation theory confirming, for the weak-coupling limit, validity of 
the second law. (Another possibility of deriving basic kinetic equations (|l^) below 
may be connected with the Tokuyama - Mori identity stemming from the Heisenberg 



equations of motion for quantum operators - see, for another model, e.g. pO|.) We 



fully rely upon the Davies mathematics and use, for the chosen model (i.e. a specific 
case) no additional approximations. Hence, rigour of our approach below is that one 
of the general Davies theory. We only deviate from a subsequent standard application 
of general Davies theorems to the weak coupling situation by just another, but equally 
admissable and physically motivated, choice of the scaling parameter. This choice 
makes the theory physically applicable also beyond the weak coupling limit. 

In order to convince sceptics concerning physical applicability of the scaling ideas, we 
also present completely independent mathematical non-scaling arguments. Those may 
be found in the Appendix and are important in particular in the relevant long-time 
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limit. ^ 



As we have numerically verified and as it is also argued below, our results well coincide with 
those of the weak-coupling theories in the overlap region with the weak-coupling regime. 
We thus have no doubts on the validity of the second law there. However, beyond the weak 
coupling regime, our results become appreciably different from standard ones. Physically, 
we have reasons supported by arguments to understand that: This is a deviation from the 
canonical state of the system caused by its non-negligible coupling to the bath. As far as 
the underlying mathematics is concerned, it is general (valid for any choice of the scaling 
parameter) and cannot be consequently sometimes correct and sometimes not. It can be 
either correct or not in general; no other alternative exists. The first alternative provides 
solid basis for arguments in favour of correctness of our approach while the second one (that 
was even never suggested or indicated) deprives even the weak coupling relaxation theory 
of its mathematical foundation. Rejecting these alternatives would 

• either mean to question physical as well as mathematical principles (including the Li- 
ouville equation) on which all the existing renown of description of kinetic phenomena 
via corresponding kinetic equations (depending on the regime in question) relies, 

• or to admit that the very principles of quantum mechanics of the open system would 
have to be complemented by, e.g., some additional requirements not admitting models 
of the type investigated here. 

In view of existing experimental evidence in favour of the quantum theory as well as because 



of experimental results indicating violations of the second law in experiment [|Tl|, 14, [2^ (see 



also for theoretical interpretation of the experimental results reported in |jT^ as well as 



previously '^^), both the latter possibilities seem, in our opinion, rather unlikely. In any 



^ Doubts recently appeared about ability of the Davies scaling theory to yield reliable long-time predictions 
on the relaxation p^ . The exactly solvable model considered was, however, bilinear in creation and 
annihilation operators. Such models, like that of harmonic phonons, are known on the other hand to 
be unable to describe diffusion regime underlying, e.g., the effect discussed in Thus, should be 
taken only as a warning against thoughtless application of formal conclusions of the Davies theory to, 
e.g., diffusion flows where no such flows exist. Here, it means another formal argument for supporting 
conclusions of the (otherwise broadly tested) Davies theory in our main text by independent arguments 
in the Appendix. 
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case, the conclusions suggest that there is at least something in physics beyond the second 
law what is at present still insufficiently understood. 



Our system, in accordance with the experimental plasma system |TT|], is assumed to 
consist of three sites, designated as 1, 2, and 3. The reader is referred to |TT|| or ||2T[] if he/she 
is interested in the motivation for construction of the model. The latter is, as compared 
to pi|, only slightly modified here so that it now better corresponds to the experimental 



system of [|r^ . Shortly, the above sites correspond to walls of the plasma container, plasma, 
and the probe. Hamiltonian of the system reads 

3 



Hs = ^ dja^aj + J{a\a2 + a^ai) + K{a\az + 0302), (1) 
i=2 



where the zero of the energy is taken to be at the walls (site 1). Though it is not in principle 
important, we assume here €3 > £2 > 0. This corresponds to the experimental situation 
|TT| . We assume only one electron in the system that is elastically transferred between 
sites 1 and 2, and simultaneously between sites 2 and 3. This is in accordance with, e.g., 
the standard theory of the Richardson-Dushman thermal emission (for the 1-2, i.e. wall - 
plasma transfers) that is based on the idea of prevailingly elastic transfer upon, e.g., electron 
leaving surface of solids. 

The load between the probe and walls in the Sheehan's experimental set-up ||TT| is the 
location at which the electron can inelastically scatter. This means, in our case, phonon- 
assisted 3^1 transitions. The phonons involved are assumed to be those of the load, here 
designated as bath II. In addition to that, we assume another bath, say bath I, formed by 
phonons (physically, those from the walls) interacting site-locally with the electron located 
on site 1. This means that Hamiltonian Hb of the bath of our model reads as 

Hb = H^B + H^B ■> 
Hi = Y^ hu;,b%, H'J = hu^BlB^. (2) 
The electron-bath coupling, Hs-b, is given by 

Hs-B = Hg_^ + Hg^_B, 
Hs-B = -tTt H h^KG^{a\ai + a\a^){h^ + hi). (3) 



(Contingent terms with ki = K2 in Hg^^ could be turned below to just a temperature- 
dependent renormalization of site-energy ei of site 1 that we set zero here.) Here the anti- 
or commutational relations between creation and annihilation operators for electrons and 
phonons are as usual. Also designates the number of phonon modes; here it is understood 
that +00. Notice that Hg_^ is quadratic in the phonon operators; this is certainly 

admissable but contrasts with both tradition and the treatment of phonon operators in 
^s-B- This non-standard assumption is employed to preserve the finite dephasing and local 
electron heating after the Davies scaling procedure (preserving formally just the lowest- 
order effects in the scaling parameter) which now follows. Otherwise, we would have to 
involve higher-order effects in treating the dephasing what could make the theory and the 
final statements, in eyes of a sceptic reader, rather ambiguous. In any case, we have a 
freedom to choose the model as above. Concerning Hg_Q, one should also notice that the 



on-site dephasing is here, in contrast to the model from |21], on site 1. This corresponds, in 



the experimentally tested system of [TT|, to electron heating inside the walls of the plasma 
container. 

One should first of all realize that the phonon-assisted transfers 3 ^ 1 as provided by 
Hg^_Q enable the electron to move in a circle 1^2^3^1,or vice versa 1 ^ 2 ^ 3 ^ 1. 
These two circular motions cannot, however, compensate each other. This is the first of 
two basic physical observations on which the present model relies. The point is that all 
the inter-site transfers involved are elastic, i.e. symmetric, except for the 3 1 one. Such 
elastic transfers lead to a tendency of equilibration of site occupation probabilities. For 
instance, assume for a while that we had only a dimer composed of sites 1 and 2, with 
the coherent (i.e. elastic) hopping term J(a|a2 + cil'^i) decoupled from any bath. The 
Hamiltonian reads then Hs = 62*^2 + J{aia2 + H.C.) as above. Stationarity of the solution 
we are interested in implies pi2 = P2i- Contingently nonzero values of these site-off-diagonal 
elements are connected with a 1-2 bonding. If we add a mechanism breaking such bonds (but, 
for simplicity, causing no additional 1^2 transfer), magnitude of pi2 = P21 would become 
suppressed below their maximum value given by positive semi-definiteness of p. Concerning 
the site-diagonal elements p^j, their stationary values can be investigated by generalized 
master equations where memory functions determine the 1^2 balance. These memory 
functions (and also their time integrals whose ratio determines pii(+C)o)/p22(+cx3)) contain 
two channels |2^], sometimes interpreted as phonon- (or bath-) assisted and quasicoherent 
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one 



pO[] . The latter channel is, in contrast to the former one, symmetric what is the reason 



of generally comparable stationary values of pii(+oo) and p22(+oo). Full equilibration 
Pii(+oo) = p22(+oo) then follows from our mathematics below upon full ignoring site 3, as 
a consequence of elastic character of the 1^2 transitions. The same applies for dimer 2 
and 3 once it is separated from the rest of the system. These facts will be useful below. 
With that, it is then easy to show that the equality of populations pn ~ P33 stemming from 
above reasoning based on the elastic character of the J- and /^-induced elastic transfers 
cannot comply with equilibration of the inelastic phonon-assisted transitions 3 1 leading 
Pii > P33 or possibly even to pn ^ P33. This fact will be useful to understand the results 
obtained below. A word of warning is, however, worth already here: What is here now being 
explained are still just heuristic arguments supported by previous investigations that explain 
our motivation; true rigorous mathematics comes only below. As it follows from general 
arguments above as well as the mathematics of the next sections, we do understand why, 
e.g., the detailed balance conditions could in our system become violated. This is because 
these conditions do not apply to uphill or downhill transfers caused by elastic mechanisms. 
We, however, do not raise the question about violations or preserving these conditions here. 
Our rigorous mathematics below avoids such statements and formulations, and leads directly 
to the required results and effects investigated here. 

For that, let us return to our model (|T|-0). Between sites 1 and 3, there is an imbal- 
ance mechanism owing to spontaneous processes allowed by Hg^_Q that prefers, because 
of assumed £3 > 0, the 3 — > 1 transitions to 1 — 3 ones. This is the imbalance (and 
the only imbalance existing in our system) that makes domination of the circular motion 
1 — >■ 2 — i> 3 — 1 over the 1 2 ^ 3 <— 1 one in fact possible. Owing to the phonon-assisted 
(prevailingly down-hill, i.e. 3^1) character of the 3 1 transfer, this implies heat transfer 
to bath II. (Each transfer act 3 — 1 is connected, because of the energy conservation law, 
with emission of a phonon quantum into bath II. Similarly for the back transfer 3^1 and 
phonon absorption. If the former transfers prevail, we get the net heat flow to bath II.) The 
question is, however, where this heat could come from. The only possible answer is that 
it is from bath I. Really, dephasing at site 1 means nothing but a continuous emission and 
absorption of phonons from bath I that can provide the necessary energy transferred by the 
electron whose energy is not sharp. Simultaneously, this dephasing can break phase relations 
between amplitude of finding the electron at site 1 and those elsewhere, i.e. it breaks the cor- 
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responding covalent-type of bonds. Without the sufficiently strong dephasing, the electron 
distribution in the system (as prescribed by, e.g., the canonical density matrix to which the 
density matrix usually tends within weak-coupling theories considering coupling to bath as 
infinitesimally weak) would really contain such and fully developed J- and i^'-induced bonds. 
Thus, it would be stiff enough in the respect that irrespective of the above imbalance, no 
electron circular motion would finally appear. One can easily verify that by calculating, e.g., 
the electron flow between any two sites. We always get (by the way, in a correspondence 
with standard physical reasoning) zero mean flow in the canonical (i.e. zeroth-order in the 
coupling to the bath) state of the system. From this point of view, lack of on-site amplitude 
dephasing (or, in other words, that of partial violations of such covalent bonds re-appearing 
immediately once the system coupling to the bath is re-introduced as a source of corrections 
to the canonical form of its density matrix) is one of the greatest deficiencies of the weak- 
coupling kinetic approaches leading to such canonical distributions. For illustration, notice 
that site off-diagonal elements of the density matrix become, in the weak-coupling kinetic 
theories, asymptotically independent of the strength of the site-local coupling to the bath 
even when this type of the coupling causes bath-induced (and on different sites uncorrelated) 
fluctuations of site energies (see, e.g., formula (28) of [0). So, even when the model does 
describe it, its weak-coupling kinetic relaxation formalism is in principle unable to describe 
the dephasing and breaking of the bonds. Hence, in our case here, we must definitely go 
beyond the weak coupling kinetic theory in especially the dephasing rate at site 1. This is 
the second basic and, perhaps, the most important observation connected with the model 
and the effect we should like to describe here. (One could also pose a question why we should 
be so keen to describe and include such potentially weak corrections to, e.g. violations of 
covalent 1 — 2 and 2 — 3 bonds here. The point is that these violations provide perhaps 
just small corrections to the canonical density matrix of the system but, simultaneously, 
these corrections are the very source of the effect we are interested in.) One should add that 
in our approach involving other than the weak coupling approach below, we also take the 
limit of the infinitesimal coupling to bath, i.e. infinitesimal dephasing. On the other hand, 
we simultaneously scale also the hopping (transfer) integrals what makes the ratio of the 
in-phasing and dephasing constant. This, in contract to the standard weak coupling scaling, 
corresponds to reality and allows the above bond breaking. (In the Sheehan experimental 
plasma system |]1T|, no such 1 — 2 or 2 — 3 covalent bonding exists). Thus we have a hope. 
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and really do obtain below the effect expected. So far, of course, all these ideas provide us 
with at most a physical background of the model, its mathematical treatment, and physics 
beyond. So, let us now have a look at how these ideas work within a rigorous theory. 



III. DAVIES SCALING AND KINETIC EQUATIONS 



The scaling procedure we use is based on Davies |2^. (For independent non-scaling 
arguments fully supporting results obtained here see the Appendix.) We, however, extend 
our treatment beyond standard weak coupling theory in that we scale not only time and 



Hs-Bi but also the transfer (overlap, hopping, or resonance) integrals J and setting ||32 

t = t'/X^, Hs^bocX, JocA^ KocX"^ (4) 

{f playing the role of a new rescaled time.) As usual, we then project off the bath and 
let A ^ 0. Such a physical regime where intersite hopping (transfer) integrals determining 
rates of bath- free transfers inside system get comparable with rates of bath-assisted processes 
(transfers) can be definitely not that of the weak coupling but rather that of the intermediate 
or, in a sense, even contingently strong coupling to the bath. Technically, though the 
mathematics used is completely that by Davies ^6], we proceed simultaneously according 



to where the relevant formulae are rewritten in a physically understandable form. 
The Davies formalism starts by writing total Hamiltonian 

H = Hs + Hb + Hs-b (5) 

in form 

H = Ho + XHi. (6) 

Here, one should add that XHi oc A but that does not exclude the possibility that XHi 
includes also higher orders in A (A^ if (^) is accepted). Those who do not like this way of 
thinking could replace conditions J oc A^, oc A^ in (^) by J oc A, oc A and proceed as 
below. The final result is the same. 



Next, introduce superoperators 

1 



Co = kHo,...], C, = l[XH^,...]^X. (7) 



Finally, be 

P... = p^®TrB(...) (8) 
10 



(with TtbP^ = 1 implying = V) the Argyres-Kelley projector (projection superoperator) 
in the Liouville space of operators that act in the Hilbert space of the system and bath. 
Then the message of Davies (see Eq. (1.19) of or Eq. (14) of [^) is 



hm sup I \p{t) - e-^iCo+{C^)+^X^K)tp^Q^ 1 1 _ 
r+oo 

A^/C... = / dxTre (e''^°"(-^£i)e-'^°"(l - V){-iCi){p'' ® ...)), 

JO 

(A)... =TrB(p^®/:i...)- (9) 

(Finite constant a is here arbitrary.) Here p{t) = TiBP^^^it) is the density matrix of the 
complex 'system + bath' with its time-development determined from the exact Liouville 
equation i-j^P^^^if) = {Cq + Ci)p^'^^{t). The assumptions used were in particular 

• that the density matrix of the system and bath p^^^{t) is initially separable, i.e. that 

p^+^(0)=p^®p(0), (10) 

and 

• that VCq = CqV. This condition can be, however, well fulfilled as far as, e.g., p^ = 
f{H^, H^^), [H^, H^J] = 0. This is in particular in our case because we are forced to 
assume 

to be able to introduce properly the initial temperatures of baths / and // separately. 

Suspicious reader with potential objections about applicability of the scaling theories to, 
in particular, the long-time relaxation phenomena is again referred to the Appendix, in 
particular its concluding remarks. Here, we should just like to add one comment: Situations 
are known when (^ cannot be used for finite A's (as in Nature). Perhaps the simplest model 
of this type is that of a single particle on a periodic chain where the system is artificially 
introduced via a few chosen sites (with the particle or without it), with the rest being the 
(formal) bath and the system-bath coupling (bringing the particle to or out of the system). 
Then matrix elements of pit) between sites of the system decay algebraically |^ while 
g-t{£o+(£i>+«A2/c)tp^g^ would indicate exponentially dominated decay. The model by Novotny 
28| is also of this type. With our real bath known to yield a well-defined relaxation, we 
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shall for the sake of brevity here postpone discussion of such singular cases with artificially 
introduced baths to another publication. 

Meaning of the mathematically exact statement in the first row of (|^) is that time devel- 
opment of p{t) as prescribed by the exact Liouville equation for the density matrix p^~^^{t) 
of the complex 'system + bath' is not discernable, in the scaling limit A ^ 0, from that one 
dictated by the kinetic equation for the density matrix p{t) of just the system 

^±p{t) = {Co + {C^)+^X'lC)p{t). (12) 

This general and exact result should now be specified according to the choice of Hq and XHi 
in®. 

Two main possibilities exist. 

• Either we accept so called weak-coupling scaling according to van Hove (and often 
automatically accepted in general situations) 

t = t'/X'^, Hs^B oc \, J = const, = const (13) 

(again with A 0) which would correspond to the choice 

Ho = Hs + Hb, XH,=Hs-b. (14) 



Then A^/C in (0) is nothing but the weak-coupling relaxation superoperator and ([12D 
reduces to the Redfield equation (before using the Redfield secular approximation) 



33| , p4| , |35| . The relaxation is then practically to the canonical state of the system. 
Or we assume (^ what means to identify 

Ho = Hs\j=K=o + Hb, XHi = Hs-b + Hs\j^o,k^o — Hs\j=k=o- (15) 
Then several things have to be realized: 

— Though the relaxation superoperator A^/C in ( 1T^ [as defined in (^] involves 
formally also higher orders in A, application of the rule Tr^ ( [aj„a„, .. .]) = 
[al^ttn, Tt^b{- ■ •)] yields that in fact just second order terms in A survive. 

— Because now Cq = ^[Hs\j=k=o+Hb, ■ ■ ■] and Hs\j=k=o is site-diagonal, the relax- 
ation is not any more (like in case of the weak-coupling choice ([l^ ) among eigen- 
states of Hs but, instead, among those (site-local eigenstates) of Hs\j=k=o- On 

12 



the other hand, in ([T^), terms Cq + reproduce \[Hs, . . .] = j\Hs\j^q^K) ■ ■ •] 
what is a free propagation among eigenstates of the full (site-off-diagonal) Hs- 
This competition between site-local and site off-local tendencies of the time de- 
velopment is what makes the dynamics much richer than in the weak-coupling 
case. 

— The fact that A^/C describes relaxation in the site-local basis is not owing to 
neglecting anything or any type of approximation. It is owing to choice of another 
regime; in our case that one in which J- and i^'-induced processes become at most 
comparable with those caused by the system interaction with the bath. 

The physical argument in favour of the form of the relaxation superoperator correspond- 
ing to the choice (|T5|) is that we are interested in the regime in which the bath-assisted 
processes inside the system are at least comparable with, or even dominating over the in- 
ternal transfer processes inside the system caused, in our case, by the J- and i^-dependent 
hopping terms in Hs (|l]). Once we realize that the weak-coupling theory presumes, in the 
sense of (^3]), the system-bath coupling to be infinitesimal, i.e. infinitely times weaker than 
all other relevant competing transfer and relaxation mechanisms, this excludes scaling (|13|), 
i.e. the choice (0), in the regime considered here. On the other hand, it allows to use (p, 
i.e. the choice (|T3p. That is why we shall below stick to this alternative. So, we use K 
from (j^) with (|T3|) for our model (j^-j^). This means the Redfield form of K, in the localized 
basis as a consequence of another (than the weak-coupling) physical regime, i.e. also cor- 
respondingly another form of identification of perturbation. So, this form of the Redfield 
tensor is definitely not consequence of any additional approximation applied to the Redfield 
form of the relaxation superoperator in the weak-coupling regime (for discussion of such an 
approximation in the weak-coupling regime see [^). After some straightforward algebra. 
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then turns in the site representation to 
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where Tf(/7) = l/{kBPi(n)) are the initial temperatures of Baths I and II; nB{P,z) is the 
Bose-Einstein phonon distribution function. and are the Golden Rule formulae for 
transfer rates 1^3 and 3 — 1. Note that F-f and F^ are different solely in that the 
latter involves a 1 + term, whereas the former has only ub- Physically, this corresponds 
to F| involving only bath-assisted stimulated up-hill transitions (absorption), whereas Fj^ 
involves both bath-assisted spontaneous and bath-assisted stimulated down-hill transitions 
(emission). Finally, 2F determines the rate of dephasing arising from local electron-energy 
fluctuations from Bath I, and also the rate of electron heating in the walls (site 1). 

A few comments are worth mentioning already here. First, notice that temperature T/ 
of bath I enters (|T6|) only via the dephasing (and simultaneously heating) rate 2F. This 
rate depends, however, also on strength and details of the electron coupling to bath I. 
Thus, moderate changes of Tj may be well compensated by those of the coupling and vice 
versa. As there are no abrupt qualitative changes expected with moderate changes of the 
coupling, only continuous changes of, e.g., the electron 1 ^ 2 ^ 3 flow are expected 
when (the initial) temperature Tj of bath I sinks below that (i.e. Tjj) of bath II. This is 
important for interpretation of the result to be obtained below. Finally, concerning (the 
initial) temperature Tjj of bath II: We shall assume here the inequality 

ksTjj ^ ea- (19) 

The opposite inequality would imply high-temperature regime in which the spontaneous 
3^1 processes would become negligible with respect to stimulated ones. So, asymptoti- 
cally, F|/Fj would turn to unity and the driving force in the circle 1 ^ 2 — 3 ^ 1 would 
disappear implying disappearance of the electron flow. This means that also the contradic- 
tion with the second law we aim at would disappear, in a full correspondence with the Bohr 
correspondence principle and the Martynov proof of validity of the second law in classical 
statistical mechanics |jl6|. On the other hand, the low-temperature limitation (|l^) is not 
severe. With, e.g., €3 ~ leV, temperatures Tjj appreciably higher than room temperatures 
are viable. 



15 



IV. ANALYTICAL SOLUTION AND STEADY-STATE HEAT-FLOW 



Let us henceforth investigate the stationary situation. Then the left hand side of ([I6D 
equals zero so that we have a homogeneous set of 9 linear algebraic equations for the sta- 
tionary values of the electron density matrix. The matrix rank is, however, only 8 since the 
sum of its first three rows is zero. Thus, the set can (and must) be complemented by the 
normalization condition 

pii + P22 + P33 = 1- (20) 

This provides us with a complete inhomogeneous set of 9 linear algebraic equations for 
9 elements of the particle density matrix. The site-diagonal matrix elements pjj give the 
probabilities of finding the electron at site j. Full algebraic solution of this set of equations is 
possible but unwieldy. That is why the set will be solved numerically below. First, however, 
we shall analytically prove that there is always, for nonzero temperatures T/ and T//, a 
positive electron flow 1 — 2 — 3 — 1 implying, for Ti < Tu, violation of the second law 
of thermodynamics in its Clausius formulation. The proof is made by logical contradiction. 
We stress from the outset that the mathematical derivation of our starting equations (jTB]), 
including scaling, involves no approximation; therefore, aside from the model assumptions, 
our treatment here is fully rigorous, in the full mathematical meaning of the word. 

Cyclic mean electron flow in the system (taken as positive in the direction 1 2 — 3 ^ 
1) can be written, on grounds of physical meaning of and Fj in (|18]), as 

J = T^ps3-T^pn. (21) 
From the first and second equations of (|16D one also has 

J=^J{P21-P12) = ^K{P32-P23). (22) 

n n 

These formulae can be also easily derived from elementary quantum mechanics. Assume 
now that no heat flows directly from bath 1 to bath II and vice versa. So, only the electron- 
mediated heat flow from I to II may appear. Since the 3^1 phonon-assisted transitions 
are associated with influx or efflux of energy (heat) to or from Bath II (proportionally to 
the magnitude of es), the total mean heat flow from Bath I to Bath II may be written as 

Q = esJ. (23) 
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This is the main quantity we are interested in. 

Let as now assume, in accord with our strategy of proof by contradiction, that there is 
no heat flow between the baths, i.e. 

2 = 0. (24) 
Since es > 0, this imphes, via (|2T|) and (p3D, that 

Tipss - r^pu = 0. (25) 



The reader could easily recognize that ( pSj) is nothing but a detailed balance condition for 
inelastic phonon-assisted 1 3 direct transitions. Because of ( P^ this also implies that 

P21 = Pl2, P32 = P23- (26) 

Now, summing the forth and fifth equations of ([TBI) with zero left hand side, we get 

= jipsi - P13) + i-2tT - ^r^)(pi2 + P21) + |(P21 - P12). (27) 

Similarly, from the sixth and seventh equation, and also from the eighth and ninth equation 
of (|16D (always with zero left hand side), we get 

= "^(P21 - P12) - 22r(pi3 + P3i) + ^(P31 - P13) + ^(P23 - P32) (28) 



and 



= ^(Pl3 - P3l) + ^^-^{P23 - P32) - ^r|(p23 + P32). (29) 



In combination with (^), Eqs. ( p^^ - |29|) give 



and 



P12 = P21 = ^-p , i-p psi, (30) 
psi = P31, (31) 



2J 

P23 = P32 = -^5>m P31. (32) 



Let us now take difference of the fourth and sixth equation in (16) in the stationary state. 



Owing to (|30|-pl|), it gives 
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Similarly, from the difference of the sixth and seventh equation and taking into account 
(0) and (ID, we obtain 







W + dr^ + + 2(2r + r, + 



9m P31. (34) 



As the expression in the square brackets is always positive, this implies that 

3m = 0, (35) 

i.e. using 



Pl3 = P31 = Pl2 = P21 = P23 = P32 = 0. (36) 

On the other hand, from (pSj) and ([35|), we get that in the stationary state 

Pll = P22- (37) 

One should realize that conditions (|37|) and pi2 = P21 (see (|36D) obtained so far fully corre- 
spond to what has been said above about no-flow equilibrium inside the dimer '1 - 2'. 
The eighth or ninth equation of (|T6|) yield, in the stationary state and with the help of 

(0), 

P22 - P33 = 0. (38) 

Together with (|20|) and (pSf), it provides an inhomogeneous set of three linear algebraic 
equations determining the site occupation probabilities (all the time provided that the no- 
fiow condition (P^ ) used above applies). The solution reads 

r, Tt 

= F ^OF ' ^22 = P33 = p . ^-p ■ (39) 

This result, on the other hand, contradicts (pT]). This is the required contradiction implying 
that (p^ ) cannot be correct. One can also ask what is the reason for the contradiction. 
Clearly, ( P7| ) would be satisfied by ( |5P| ) if there were Fj — = 0. That would, however, 
mean to disregard the spontaneous processes that are responsible for the difference on the 
left hand side. The spontaneous processes are, however, purely quantum. Similarly, one 
can easily observe that (^) becomes fully compatible with (|37D in the limit of the infinite 
temperature T// — * +00. The infinite temperature limit means, however, the classical limit 
(the Bohr correspondence principle). All that is why we can understand the violation of the 
second law we arrive at below (as well as in other models yielding such a striking conclusion 
- see above) as a consequence of quantum effects. 
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So, there is always an electron circular flow in the system implying (not in general but) 
in our speciflc situation nonzero heat transfer Q (as given by ( p3D and (^) or (0)) between 
baths I and II. The last questions to be solved before we resort to a numerical study are 
what is its orientation and how the conclusion contradicts the second law. 

V. VIOLATION OF THE SECOND LAW 

In order to infer what is the orientation of the mean heat flow, let us turn to above 
formulae ( pi]) and (pSf). From (|23D we get that signs of Q and J' coincide (63 > 0). As for 
the latter, we remind that P33 is always (as a site occupation probability) positive and that 

disappears for Tjj 0. Thus, from (0), we get that is, in the low temperature limit 
of bath II but arbitrary nonzero Tj, always positive. This is, by the way, also what our 
numerical results show. 

Let us now increase T//. One should realize that Q is a continuous function of Tjj and 
never turns to zero. (For that, see the above proof.) So it should remain positive even when 
Tjj becomes greater than Tj. (In fact, owing to intermixture of Tj with details of coupling 
to bath I in Hg_Q inside F, nothing can happen at the moment when T// passes T/. This 
fact was also conflrmed numerically.) Positive values of Q mean, however, a positive rate 
of heat transfer from bath I to bath II which thus goes, for T/ < T//, against temperature 
step. As the heat transfer is spontaneous (there is no external expenditure of energy or 
whatever else conditioning this transfer), this conclusion explicitly contradicts the Clausius 
form of the second law stating that such processes are impossible. On the other hand, 
the conclusion obtained analytically here (and verifled numerically below) that the second 
law is in our system really violated fully corresponds to conclusions of |]1T| where, for an 
experimental plasma system corresponding to the above model, the universal validity of the 
second law was flrst seriously challenged. 

VI. NUMERICAL RESULTS 

In order to verify the above conclusions, we have solved the set ([I6D and (pO|) numerically. 
There was also a secondary reason for this numerical study: Analytically, we were unable 
to prove that the heat transfer really turns to zero in the limit of zero temperature T/ of 
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bath I. This is what must be expected physically because in such a limit, there is no heat 
available in bath I to be transferred to bath II. In just other words: There is no dephasing 
in this limit between sites 1 and 2. So, the covalent bond 1-2 should become perfect, making 
thus the electron (and consequently also the heat) flow impossible. 

Fig. 1 shows typical results. We designate 70 = r-|-[e^^^^^ — 1]. Three things are worth 
noticing: 

• In accordance with the above analytical arguments, the mean heat flow Q is always 
positive (i.e. going from bath I to bath II). 

• With decreasing dephasing rate 2T corresponding to decreasing temperature Tj, Q 
turns apparently to zero. 

• For constant rate 2r, Q is only very little dependent on temperature Tjj of bath II. 
Slight increase as well as decrease with T// are both possible. This may be interpreted 

result of two competing tendencies: 

— Increasing Tjj increases also the rate of dephasing between sites 2 and 3 caused 
by nonzero and T/j-dependent terms —^^1 3-3 and 4-4 elements of block C 
in (|16|). Similarly the terms — in 4-4 and 5-5 elements of block A in ( p!6[ ) 
contributing to dephasing of sites 1 and 2. This leads to greater violations of the 
2-3 and 1-2 bonds, i.e. to increase of Q. 

— Increasing Tjj on the other hand implies relatively decreasing role of the sponta- 
neous processes in 3 1 transitions what means suppression of Q. 

Dependence of Q on Tu is, however, always very small. 

VII. CONCLUSION 

We have obtained a spontaneous heat flow between two macroscopic baths that is owing 
to a speciflc activity of our microscopic single-electron system not aided from outside. So, 
starting from rigorous mathematics of the quantum theory of open systems, a contradic- 
tion with the second law of thermodynamics has been obtained for the model in question. 
Remind that except for model assumptions, no approximations were made that could be 
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made responsible for the effect, and that the model corresponds to an experimental system 



positively tested in [11]. This indicates that one should choose between just two alternatives: 



There is still something hidden in physics beyond the second law what is at present 
not fully understood. This possibility might also mean complementing contemporary 
quantum mechanics and present philosophy of quantum-mechanical modelling in order 
to reconcile the quantum theory with (presumably) universally valid thermodynamics. 

The second alternative is to refrain from the so far universally assumed validity of 
thermodynamics in the macroworld. One should realize that though our system is 
microscopic, appending the macroscopic reservoirs turn the physics to the macroscopic 
one. 
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IX. APPENDIX: DERIVATION OF THE STARTING EQUATIONS BY NON- 
SCALING ARGUMENTS 

Arguments in favour of violation of the second law as reported above critically depend 
on the existence of the heat flow from bath I to bath II, i.e. on the existence of the electron 
mean circular flow J . Existence of the flow has been above proved using scaling which 
does not belong to a generally accepted weaponry in kinetic theories. Moreover, though 
standard practice conflrms such a possibility, one could ask about justiflcation to use such a 
(as well as any other) kinetic approach beyond (time) limits of the kinetic regime. (Notice, 
e.g., that constant a in ... supo<A2t<a ••• in (|^) is always flnite.) That is why we present 
another treatment below that is non-scaling but fully conflrms the above conclusions. This 
treatment is then fully resistive even against such objections as it assumes, as experiments 
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require, taking first (though after the thermodynamic hmit of the bath) the dc, i.e. the 
infinite time limit. Only then (if at all) discussion based on smallness of individual terms in 
the Hamiltonian (coupling constants etc.) comes into question. 

For reasons connected with unreliability of the finite-order approximations in convolution 
theories discussed in we refrain from convolution theories of the Nakajima-Zwanzig 
PH , type. Instead, we adopt the formalism based on the time-convolutionless (i.e. time- 
local) Generalized master equations (TCL-GME). These were first suggested by Fulihski and 



Kramarczyk [^, 44 but more operative (and in fact equivalent [|^) are those by Shibata, 
Hashitsume, Takahashi, and Shingu [|^, The starting Shibata, Hashitsume, Takahashi, 
and Shingu identity (derived as a direct consequence of the Liouville equation) reads 



Vp^+^{t) = -iVC[l + i f exp{-z(l - V)Ct}{1 - V)CVexp{iCT} dr]'^ 

Jo 

■ [exp{-z(l - V)Ct}{l - P)p^+^(0) + Vp^+^'it)]. (40) 



Now we make two steps. First we take for V the Argyres-Kelley projector (P), and assume 
the initial condition (10). Let us now split our Hamiltonian as in (^) and introduce Cq and 



XCi as in (0). This reduces, for both the weak coupling ([l^) and our identification (|15|), 
equation (|40|) to 

|p(t) = {-^C, - ^{C,) + A^/C^^^~^*^^(t))p(t). (41) 
Here, because VCq = CqV, 
y2^TCL-GME^^^ . . . = -Ttb {Ci[l + t j\xp{-t{l - V)Ct}{1 - V)CiV exp{tCT} dr]-^ 

■ exp{-i(l - V)Ct}{1 - V)CiVexp{iCT} dr{p^ ® . . .)) . (42) 

Notice that the relaxation superoperator X^lC^'-^^~'-^'^^^ (t) still involves also higher-than- 
second order terms in A. It is also time-dependent. This time dependence is determined 
by the decay-to-zero of the integrand (after the implicit but everywhere assumed thermo- 
dynamic limit of the bath). The characteristic times are given by dephasing in the bath. 
Once we disregard a transient (and for our purposes fully unimportant) initial time period, 
we may turn time t in (^) to infinity. 
We should like to stress that 

• turning t to infinity in A^/C^'"^~*^*^'^()f:) still does not mean treating just the final result 
of the finished process of relaxation; 
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• the long-time limit taken in p^) is dictated by the experimental situation (we are not 
interested in short-time transient relaxation effects) and time t could be, beyond some 
bath dephasing time, well taken arbitrarily large, in accordance with the experiment. 

• This is unlike the situation with the value of A. In Nature, this coupling constant is 
determined, for a given experiment, once for ever and cannot be arbitrarily changed. 

That is why one should expand, if at all, in powers of A only after the long time limit 
t +00 in \^lC^'^^~'^^^^{t). That is what we do here. If we take t +00 in (H) and then 
limit our attention to the lowest nonvanishing (i.e. the second) order in A, (|4^) then turns 
to 

= -Ttb (^Ci I^"^ exp{-iCor}{l - V)CiV exp{iCor} drip"" . . + 0{\^). (43) 

The first term on the right hand side of (|4^) should now be compared with that of A^/C in 
d). The point is that if identity of these two expressions for the relaxation superoperator is 
established, the TCL-GME (|iT| ) becomes, for both the standard weak coupling choice (|T3|) 
and that one by Capek and Barvfk (p!5[), fully equivalent to ([T^). 



At the first sight, there is a similarity but no identity observed. In order to discuss 
the point in detail, let us make several physically motivated and justifiable steps. First, the 
integrations, in both A^/C in (^ and \^K7'~^^~^^'^^{t) in (^), lead (upon explicit introduction 
of the matrix elements involved) to distributions lim5^o+ /^E^-ins ^^^^ approximate as 
7Th6{AE) {AE being relevant differences of eigenenergies of Hq). This step amounts to 
neglect term v.p.^ (here v. p. means the fraction in the Cauchy sense); such imaginary 
terms are standardly interpreted as just renormalizations of transfer (hopping or resonance) 
integrals. (In this connection, notice the imaginary unit i in the {—iCo —i{Ci)) term on the 
right hand side of (^T]) or (P).) Doing so, the mn- matrix element of A^/C ... in (^ turns to 



(A^/C . . .) =-l J2T.T.{M^Hi\sa){sa\XHi\qK){qi^\p^ ® . . . |n/i)7r5(E,, - Em^) 

\ / mn n 



— {mii\\Hi\sa) {sa\p^ ® . . . \qK){qK\\Hi\np)7rS{Esa — Em^) 

— {mp\XHi\qK) {qK\p^ ® . . . \sa){sa\\Hi\np)Tc5{En^ — Ego) 
+ {mp\p^ ® . . . \qK,) {qK\XHi\sa) {sa\XHi\nfj,)TT5{Enfr — Ega)} 
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+ (X^A)C..) . (44) 

V / mn 



Here, we have used the notation \mfi) = |m) \fi) etc. where \m) and \fi) are respectively 
eigenstates of Hq — Hb and Hb', = + are the corresponding eigenenergies. In 
the same way, disregarding aheady the transient time- dependence and higher-than-second 
order terms {by the definition absent in @) 



^2j^TCL-GME^ \ ^ -lY.Y.T.{i^^^\>^H^\sa){sa\\H^\qK){qK\p''®. . . \nf^)7r5{E,^-E, 

/ mn h '-^ '—' '—^ 



qn ) 



(m/i|A-f/'i|scr)(scr|p^ ® . . . \qK) {qK\\Hi\nn)Tr5{Egf^ — E. 



nfij 



-{mfi\\Hi\qK) {qK\p^ ® ... \sa) {s(7\\Hi\nfi)'JT6{Em^, - Eg^) 
+ {mp\p^ . . . \qK,) {qK\XHi\sa) {sa\XHi\nfi)TT6{Eqf, — Eg^)} 

+ (\'AIC^^'^-^''^...) . (45) 

V / mn 

Obviously, the second and third terms from (|^) equal respectively to the third and second 
terms in (^). As for the remaining terms, there are differences between the Davies ( ^4|) and 
TCL-GME (^51) results that may become important when, e.g., different matrix elements of 
the coupling to phonons interfere. 

In order to show that this is not our case for the above model and the regime investigated, 
let us make several specifications and observations. 

• First, assume that 

(/i|p^|z/) = 5^,,p^. (46) 
This is consistent with assumptions od both the scrutinized approaches. 

• Realize that for the model specified by (^^, the two terms in the system-bath coupling 
then do not interfere. For the site-local coupling, it is trivial to see the equivalence of 
the first and the fourth terms in (|4^ ) and (|45|). So, we turn our attention to just the 
site off-local coupling term Hg^_B in (^. 

• We specify our reasoning here to the regime corresponding to our choice (|1^). Then 
Hs\j=K=o is site-diagonal, i.e. the Latin summation indices in (HI) and (H^) are sites. 



Because of the last point, equivalence between the first rows of (5^) and ( |^ ) and, similarly, 



between the fourth rows of ( PD and (^3[), is then for the site off- local coupling easily seen. 
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and 



[For example, with XHi being in the first row of (PD and (|45|) substituted by Hg^_^ from 
(0), the form of Hg^_Q imphes that nonzero contribution appears just for m = q. Because of 
(^), only terms with k, = ^ contribute what makes the equivalence between the first rows 
of ( PD and explicit.] As for the last terms in ( PD and (^^) ignored so far, we get 

( A'A/C . . . ) = / dxTvB (e'^^^'rie-'^^'P^ifp^ ® ...)) 

V /mn Jq 

= j-Yl '^{{mfi\\Hi\aa) {a\p^\fi) {au\[\Hi, (g) . . .]\nu)']TS{Eaa - Em^) 

- {mu\[\Hi,p^ (g) . . .]\au) {p\p^\a) {aa\\Hi\np)7iS{Enf, - Eaa)}- (47) 

1 c+co 

= -j dxJ2 J2[{mp\XHi\aa){a\e-'''^'^/''p''e'''^''/^\p) 

" fiuafS ah 

^^^^^^ -i{aiy\XHi\bf3){bf3\p^ . . . \nu) ^ i{au\p^ . . .\b(3){bf3\XHi\nu) 

Eau — Ebi3 — iO+ Ebp — Eny — iQ + 

^-i{mv\\Hi\b(3){b(3\p^ ® ... \av) ^ i{mv\p^ ® ... \bl3){bl3\\Hi\av) ^ 

Emu — Ehfs — iO+ Eh0 — Eav — «0+ 

• {ii\e-'^^''l^p^e'^^''/^\a){aa\\Hi\np)] + C(A^). (48) 

Clearly, expressions (^) and ( ^S]) are discernibly different even when we omit the renormal- 
ization terms. Fortunately, with ( ^61) and the fact that both the terms in Hs-b in (0) are 
off-diagonal in the phonon indices, both ( P7| ) and ( ^SD are in fact exactly zero. This makes 
the proof of full equivalency, for our model and to the lowest perturbational order, of the 
relaxation superoperator K, as derived from the scaling Davies theory with that one derived 
by the non-scaling time-convolutionless Generalized Master Equation theory complete. The 
latter theory thus provides independent non-scaling way to our above equation ([T6| ) with 
(0) forming basis of the above discussion. The important point is that we have here, af- 
ter taking the thermodynamic limit of the bath, first turned the real physical time behind 
dephasing time of the bath, i.e. potentially even to infinity. Only then we have discussed 
the form of the relaxation tensor determining the relaxation process up to infinite times as 
valid for small couplings. Hence, unlike the scaling theories, strength of the coupling never 
comes into any competition with time limitations of the theory so that no objections can be 
raised that the kinetic theory used is inapplicable behind some long critical times. 
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Figure captions 



Figure 1: Spontaneous energy flow Q in units 4K'^/h from batli I to batfi II as a function 
of temperature Tu = l/{kBPn) of bath II. We set J = X = 0.5 eV, ea/X = 4, eg/X = 2, 
h'fo/K — 0.02, and hV/K — 10~^, 10~^ and 10~^ (decreasing dephasing rate, i.e. also 
heating, at site 1) for curves a), b), and c), respectively. Notice that 2r incorporates also 
temperature T/ of bath I. 
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